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We consider a way to generate operational inequalities to test nonclassicality (or quantumness) of multimode 
bosonic fields (or multiparty bosonic systems) that unifies the derivation of many known inequalities and allows 
to propose new ones. The nonclassicality criteria are based on Vogel's criterion corresponding to analyzing the 
positivity of multimode P functions or, equivalently, the positivity of matrices of expectation values of, e.g., 
creation and annihilation operators. We analyze not only monomials, but also polynomial functions of such 
moments, which can sometimes enable simpler derivations of physically relevant inequalities. As an example, 
we derive various classical inequalities which can be violated only by nonclassical fields. In particular, we show 
how the criteria introduced here easily reduce to the well-known inequalities describing: (a) multimode quadra- 
ture squeezing and its generalizations including sum, difference and principal squeezing, (b) two-mode one-time 
photon-number correlations including sub-Poisson photon-number correlations and effects corresponding to vi- 
olations of the Cauchy-Schwarz and Muirhead inequalities, (c) two-time single-mode photon-number correla- 
tions including photon antibunching and hyperbunching, and (d) two- and three-mode quantum entanglement. 
Other simple inequalities for testing nonclassicality are also proposed. We have found some general relations 
between the nonclassicality and entanglement criteria, in particular, those resulting from the Cauchy-Schwarz 
inequality. It is shown that some known entanglement inequalities can be derived as nonclassicality inequal- 
ities within our formalism, while some other known entanglement inequalities can be seen as sums of more 
than one inequality derived from the nonclassicality criterion. This approach enables a deeper analysis of the 
entanglement for a given nonclassicality. 

PACS numbers: 42.50.Ar, 42.50.Xa, 03.67 .Mn 



I. INTRODUCTION 

Testing whether a given state of a system cannot be de- 
scribed within a classical theory, has been one of the fun- 
damental problems of quantum theory from its beginnings 
to current studies in, e.g., quantum optics fil 4l2l|. c ondensed 
matter (see, e.g., Refs. fl [l3ll ). nanomechanics Ill4l[l5ll . and 
quantum biology (see, e.g., Ref. Oil). Macroscopic quantum 
superpositions (being at the heart of the Schrtidinger-cat para- 
dox) and related entangled states (which are at the core of the 
Einstein-Podolsky-Rosen paradox and Bell's theorem) are fa- 
mous examples of nonclassical states which are not only phys- 
ical curiosities but now fundamental resources for quantum- 
information processing iTTvtl . 

All states (or phenomena) are quantum, i.e., nonclassical. 
Thus, it is quite arbitrary to call some states "classical". Nev- 
ertheless, some states are closer to their classical approxima- 
tion than other states. The most classical pure states of the 
harmonic oscillator are coherent states. Thus, usually, they 
are considered classical, while all other pure states of the har- 
monic oscillator are deemed nonclassical. The nonclassicality 
criterion for mixed states is more complicated and it is based 
on the Glauber-Sudarshan P function 111 01 ■ A commonly ac- 
cepted formal criterion which enables to distinguish nonclas- 
sical from classical states reads as follows (Lsj-lg]: A quantum 
state is nonclassical if its Glauber-Sudarshan P function can- 
not be interpreted as a true probability density. Note that, ac- 



cording to this definition, any entangled state is nonclassical, 
but not every separable state is classical. 

Various operational criteria of nonclassicality (or quantum- 
ness) of single-mode fields were proposed (see, e.g., Refs. 
R flslfisjll and references therein). In particular, Agarwal and 
Tara JM, Shchukin, Richter and Vogel (SRV) UM HI pro- 
posed nonclassicality criteria based on matrices of moments 
of annihilation and creation operators for single-mode fields. 
Moreover, an efficient method for measuring such moments 
was also developed by Shchukin and Vogel Il23ll . 

It is not always sufficient to analyze a single-mode field, 
i.e., an elementary excitation of a normal mode of the field 
confined to a one-dimensional cavity. To describe the genera- 
tion or interaction of two or more bosonic fields, the standard 
analysis of single-system nonclassicality should be general- 
ized to the two- and multi-system (multimode) case. Simple 
examples of such bosonic fields are multimode number states, 
multimode coherent and squeezed light, or fields generated 
in multi-wave mixing, multimode scattering, or multi-photon 
resonance. 

Here, we study in greater detail and modify an operational 
criterion of nonclassicality for multimode radiation fields of 
Vogel l24ll . which is ageneralized version of the SRV non- 
classicality criterion ll2ll I22I1 for single-mode fields. It not 
only describes the multimode fields, but can also be applied in 
the analysis of the dynamics of radiation sources. This could 
be important for the study of, e.g., time-dependent correlation 
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functions, which are related to time-dependent field commu- 
tation rules (see, e.g., subsections 2.7 and 2.8 in Ref. J3l). 

A variety of multimode nonclassicality inequalities has 
been proposed in quantum optics (see, e.g., textbooks JD-f^], 
reviews f7 Hllll . and Refs. Il25l - [4ll0 and tested experimentally 
(see, e.g., Refs. Ji^-SMl)- The nonclassicality criterion de- 
scribed here enables a simple derivation of them. Moreover, 
it offers an effective way to derive new inequalities, which 
might be useful in testing the nonclassicality of specific states 
generated in experiments. 

It is worth noting that we are analyzing nonclassicality cri- 
teria but not a degree of nonclassicality. We admit that the 
latter problem is experimentally important and a few "mea- 
sures" of nonclassicality have been proposed ll49l - l58ll . 

Analogously to the SRV nonclassicality criteria, Shchukin 
and Vogel 115 911 proposed an entanglement criterion based on 
the matrices of moments and partial transposition. This crite- 
rion was later amended ll60ll and generalized l6lll to replace 
partial transposition by nondecomposable positive maps and 
contraction maps (e.g., realignment). A similar approach for 
entanglement verification, based on the construction of matri- 
ces of expectation values, was also investigated in Refs. J6^- 

m. 

Here we analyze relations between classical inequalities de- 
rived from the two- and three-mode nonclassicality criteria 
and the above-mentioned entanglement criterion. 

The article is organized as follows: In Sect. [II] a nonclas- 
sicality criterion for multimode bosonic fields is formulated. 
We apply the criterion to rederive known and a few apparently 
new nonclassicality inequalities. In subsection lHI A| we sum- 
marize the Shchukin- Vogel entanglement criterion lD9l,l60ll . In 
subsection lHI CI we apply it to show that some known entan- 
glement inequalities (including those of Duan et al. l66ll and 
Hillery and Zubairy 16711) exactly correspond to unique non- 
classicality inequalities. In subsection lHI Dl we analyze such 
entanglement inequalities (including Simon's criterion 1I68I1 ) 
that are represented apparently not by a single inequality but 
by sums of inequalities derived from the nonclassicality crite- 
rion. Moreover, other entanglement inequalities are derived in 
subsection HHP 21 The discussed nonclassicality and entan- 
glement criteria are summarized in Tables I and II. We con- 
clude in Sect. [IV] 



II. NONCLASSICALITY CRITERIA FOR MULTIMODE 
FIELDS 

An A/-mode bosonic state p can be completely described 
by the Glauber-Sudarshan P function defined by QH2]: 

p = J d 2 a P(a,et*)\a)(a\, (1) 

where \a) = IIm=i \ a m) an d \ a m) is the rath-mode coherent 
state, i.e., the eigenstate of the mth-mode annihilation oper- 
ator a m , a denotes complex multivariable (cci, «2, &m)> 
and d 2 a = J\ d 2 a m . The density matrix p can be sup- 
ported on the tensor product of either infinite-dimensional or 
finite-dimensional Hilbert spaces. For the sake of simplicity, 



we assume the number M of modes to be finite. But there is 
no problem to generalize our results for an infinite number of 
modes. 

A criterion of nonclassicality is usually formulated as fol- 
lows |7(ill : 

Criterion 1 A multimode bosonic state p is considered to be 
nonclassical if its Glauber-Sudarshan P function cannot be 
interpreted as a classical probability density, i.e., it is nonpos- 
itive or more singular than Dirac 's delta function. Conversely, 
a state is called classical if it is described by a P function be- 
ing a classical probability density. 

It is worth noting that Criterion 1 (and the following crite- 
ria) does not have a fundamental indisputable validity, and it 
was the subject of criticism by, e.g., Wiinsche ifTlll . who made 
the following two observations, (i) In the vicinity of any clas- 
sical state there are nonclassical states, as can be illustrated 
by analyzing modified thermal states. So, arbitrarily close to 
any classical state there is a nonclassical state giving, to arbi- 
trary precision, exactly the same outcomes as for the classi- 
cal state in any measurement. Note that analogous problems 
can be raised for entanglement criteria 16 ill for continuous- 
variable systems, as in the vicinity of any separable state there 
are entangled states. 1 84] (ii) There are intermediate quasiclas- 
sical (or unorthodox classical) states, which cannot be clearly 
classified as classical or nonclassical according to Criterion 1 . 
This can be illustrated by analyzing the squeezing of thermal 
states, which does not lead immediately from classical to non- 
classical states. 

Due to the singularity of the P function, Criterion 1 is 
not operationally useful as it is extremely difficult (although 
sometimes possible ll72tl ) to directly reconstruct the P func- 
tion from experimental data. 

Recently, Shchukin, Richter and Vogel I21I l22tl proposed a 
hierarchy of operational criteria of nonclassicality of single- 
mode bosonic states. This approach is based on the normally 
ordered moments of, e.g., annihilation and creation operators 
or position and momentum operators. An infinite set of these 
criteria (by inclusion of the correction analogous to that given 
in Ref. 16011 ) corresponds to a single-mode version of Crite- 
rion 1. 

Let us consider a (possibly infinite) countable set F = 
(A, A, ...,fi,...) of M-mode operators = /;(a, a t ), 
each a function of annihilation, a = (61, 0,2, clm), and 
creation, aJ, operators. For example, we may choose such 
operators as monomials 

M 

fi= W^ir—^, (2) 

m— 1 

where i stands in this case for the multi-index i = 
(«i,i2, •••,*2m); but the /j's can be more complicated func- 
tions, for example polynomials in the creation and annihila- 
tion operators. 
If 

/ = (3) 
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where c, are arbitrary complex numbers, then with the help of 
the P function one can directly calculate the normally ordered 
(denote d by : :) mean values of the Hermitian operator pf as 
follows El El: 



Pf:) = / d 2 c* \f(oc,a*)\ 2 P(cx,cx*). 



(4) 



The crucial observation of SRV 112 111 in the derivation of their 
criterion is the following: 

Observation 1 If the P function for a given state is a classical 
probability density, then (: pf :) > for any function f. 
Conversely, if (: p f :) < for some f, then the P function is 
not a classical probability density. 

The condition based on nonpositivity of the P function is usu- 
ally considered a necessary and sufficient condition of non- 
classicality. In fact, as shown by Sperling IT74ll . if the P func- 
tion is more singular than Dirac's (5-function [e.g., given by 
the nth derivative of 8(a) for n = 1,2, ...], then it is also 
nonpositive. 

With the help of Eq. (01, (: pf :) can be given by 



(■■Pf--) = 



in terms of the normally ordered correlation functions 

M^(p) = Tr : p), 



(5) 



(6) 



where the superscript (n) (similarly to : :) denotes the normal 
order of field operators. In the special case of two modes, 
analyzed in detail in the next sections, and with the choice of 
fi given by Eq. (|2), Eq. (O can be simply written as 

M^f(p) = Tr [ : (a^ 1 aH*H u )\a^ aH^V*) : p], 

( 7 ) 

where a = ai and b = a.2- It is worth noting that there is 
an efficient optical scheme i23ll for measuring the correlation 
functions ©. 

With a set F = (/i, fa, . . , , fi, . , .) fixed, the correlations 
© form a (possibly infinite) Hermitian matrix 



M^{p) = [M^\p)]. 



(8) 



In order to emphasize the dependence of ([8]l on the choice of 
F, we may write M^\p). Moreover, let [M^ (p)] r , with 

F 

r = (fx, ... , rjy), denote the N x N principal submatrix of 
M^ n > (p) obtained by deleting all rows and columns except the 
ones labeled by ri , . . . , risr ■ 

Analogously to Vogel's approach Q24ll . by applying 
Sylvester's criterion (see, e.g., ll60l l75ll ) to the matrix (0, 
a generalization of the single-mode SRV criterion for multi- 
mode fields can be formulated as follows: 

Criterion 2 For any choice of F = (/i, fi, . . . , f, . . .), a 
multimode state p is nonclassical if there exists a negative 
principal minor, i.e., det[M<?\p)] r < 0, for some r = 
(ri, . . . , r N ), with 1 < ri < r 2 < ■ ■ ■ < r N . 



According to Vogel 02411 . this criterion (and the following Cri- 
terion 3) can also be applied to describe the nonclassicality of 
space-time correlations and the dynamics of radiation sources 
by applying the generalized P function: 



P(a,a*) = ( oY[S(a l -a t y c 



(9) 



where a = (a\, ...,olm), with ai = ai(ri,ti) depending on 
the space-time arguments (r^, i,). By contrast to the standard 
definition of P function, symbol ° ° describes both the normal 
order of field operators and also time order, i.e., time argu- 
ments increase to the right (left) in products of creation (an- 
nihilation) operators |4J]. As an example, we will apply this 
generalized criterion to show the nonclassicality of photon an- 
tibunching and hyperbunching effects in Appendix C. 

Note that Criterion 2, even for the choice of /,; given by 
Eq. (f2]i and in the special case of single-mode fields, does 
not exactly reduce to the SRV criterion as it appeared in 
Ref. ll22ll . To show this, let us denote by M^\p) the sub- 
matrix corresponding to the first N rows and columns of 
M( n > (p). According to the original SRV criterion (Theorem 3 
in Ref. ||22ll ), a single-mode state is nonclassical if there exists 
an N, such that the leading principal minor is negative, i.e. 
det Mpj(p) < 0. Such formulated criterion fails for singular 

(i.e., det M^f 1 (p) = 0) matrices of moments, as explained in 
detail in the context of quantum entanglement in Ref. l60ll . 
Considering [M^\p)] r is equivalent to considering the 

r 

correlation matrix corresponding to a subset F' C F, with 



..,f rN ), i.e., [Mfip)} 



M$>(p). We 



F' = C/W,/, 

note that the subset symbol is used for brevity although it is 
not very precise, as the i^s are ordered collections of opera 
tors. 

Thus, by denoting 



M^{p)^[M^{p)]t 



(■■ flX :) (: ftjr 2 :) 



('■ frifjN ') \ 
('■ fr 2 f r N ') 



,(10) 



V ( : fr N fri '■) ('■ fr N fr-2 '■)''' ('■ fr N fr N ■) I 

and its determinant 

det Mf{p), (11) 
we can equivalently rewrite Criterion 2 as: 

Criterion 3 A multimode bosonic state p is nonclassical if 
there exists F, such that (p) is negative. 

F 



This can be written more compactly as: 



p is classical => \/F 
p is nonclassical <^ 3F 



df(p)>0, 



d ( ~'(p)<0. 



(12) 
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In the following, we use the symbol < to emphasize that a 
given inequality can be satisfied only for nonclassical states 

cl 

and the symbol > to indicate that an inequality must be sat- 
isfied for all classical states. 

Let us comment further on the relation between Criteria 2 
and 3 and the SRV criterion (in its amended version that takes 
into account the issue of singular matrices). Criterion 3 cor- 

(n) 

responds to checking the positivity of an infinite matrix 

defined as in © with the fa's chosen to be all possible mono- 
mials given by Eq. ©. Considering the positivity of larger 
and larger submatrices of this matrix leads to a hierarchy of 

(n) 

criteria: testing the positivity of some submatrix M N leads 
to a stronger criterion than testing the positivity of a subma- 
trix Mj^) , with N' < N. Nonetheless, when one invokes 
Sylvester's criterion in order to transform the test of positivity 
of a matrix into the test of positivity of its many principal mi- 
nors, it is arguably difficult to speak of a "hierarchy". Indeed, 
because of the issue of the possible singularity of the matrix 
we cannot simply consider, e.g., leading principal minors in- 
volving larger and larger submatrices. 

As regards the general formalism, of course by adding op- 
erators to the set F, and therefore increasing the dimension of 

in) 

the matrix Mi , one obtains a hierarchy of matrix conditions 
on classicality. Nonetheless, also in our case when moving to 
scalar inequalities by considering determinants, we face the 
issue of the possible singularity of matrices. Motivated also 
by this difficulty, in the present article we do not focus so 
much on the idea a hierarchy of criteria, but rather explore 
the approach that by using matrices of expectations values it 
is possible to easily obtain criteria of nonclassicality and en- 
tanglement in the form of inequalities. As already explained, 
this is done by referring to Observation[T]and considering fa's 
possibly more general than monomials, e.g., polynomials. 

Indeed, when we choose a set of operators F = 
(A; /2, • ■ • )' we compute the corresponding matrix of expec- 
tation values, and we check its positivity, what we are doing is 
equivalent to checking positivity of, e.g., (: f* f :) for all /'s 
that can be written as a linear combination of the operators in 
F: f = J2i c ifi- As polynomials can be expanded into mono- 

(n) 

mials, it is clear that checking the positivity of a matrix M . ' 

F 

with F consisting of polynomials, cannot give a stronger cri- 

(n) 

terion than checking the positivity of a matrix AP~ ', where 

F 

F' is given by all the monomials appearing in the elements 
of F. Of course, to have a stronger matrix criterion of classi- 
cality we pay a price in terms of the dimension of the matrix 
which is larger than Further, as we will see, 

F F 

by considering general sets F — that is, not only containing 
monomials — one can straightforwardly obtain interesting and 
physically relevant inequalities, which may be difficult to pin- 
point when considering monomials as "building blocks". It is 
worth noting that the possibility of using polynomial functions 
of moments was also discussed in Ref. |59[| in the context of 
entanglement criterion. 

Finally, we remark that to make the above criteria sensitive 



in detecting nonclassicality, the fa must be chosen such that 
the normal ordering is important in giving M^ n > . In particular, 
assuming this special structure for the fa's, there must be some 
combination of creation and annihilation operators. On the 
contrary, the inclusion of only creation or only annihilation 
operators would give a matrix M^ n > positive for every state, 
thus completely useless for detecting nonclassicality. 



A. Nonclassicality and the Cauchy-Schwarz inequality 



The Cauchy-Schwarz inequalityJCSI) for operators can be 
written as follows (see, e.g., Ref. ||5|]): 



(A*A)(&&) > \{A^B)\ 



(13) 



where A and B are arbitrary operators for which the above 
expectations exist. Indeed, (A^B) = Tr (pA^B) is a valid 
inner product because of the positivity of p. Similarly, one 
can define a valid scalar product for a positive P function. In 
detail, by identifying A = /i(a, a^) and B = /2(a, a^), one 
can define the scalar product 



<: flfj ■) = / d 2 « /*(a,a*)/,(a,a*)P(a,a* 



Then, a CSI can be written as: 



(■■ flfi ■■)(■■ fU2:)>\(: f I h 



(14) 



(15) 



Such CSI, for a given choice of operators fa and fa, can be 
violated by some nonclassical fields described by a P func- 
tion which is not positive everywhere, that is such that ( TT4| ) 
does not actually define a scalar product. We then say that 
the state of the fields violates the CSI. The nonclassicality of 
states violating the CSI can be shown by analyzing Criterion 3 
for F = (fa, fa), which results in 



/7/i:> (■■flh--) 

hfl) (■■flh-.) 



n ^o. 



(16) 



B. A zoo of nonclassical phenomena 

In Table I, we present a variety of multimode nonclassical- 
ity criteria, which can be derived by applying Criterion 3 as 
shown in this subsection and in greater detail in Appendices 
A-C. 

In the following, we give a few simple examples of other 
classical inequalities, which — to our knowledge — have not 
been discussed in the literature. In particular, we analyze in- 
equalities based on determinants of the following form: 



D(x 7 y,z) 



1 XX* 

x* z y* 
V * 



.v 



(17) 
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TABLE I: Criteria for single-time nonclassical effects in two-mode (TM) and multimode (MM) fields, and two-time nonclassical effects in 
single-mode (SM) fields. 



Nonclassical effect 



Criterion 



Equations 



MM quadrature squeezing 
TM principal squeezing of Luks et al. f321 
TM sum squeezing of Hillery f331 
MM sum squeezing of An-Tinh f391 
TM difference squeezing of Hillery f331 
MM difference squeezing of An-Tinh f401 
TM sub-Poisson photon-number correlations 
Cauchy-Schwarz inequality violation 

TM Cauchy-Schwarz inequality violation via Agarwal's test f3 11 
TM Muirhead inequality violation via Lee's test f341 
SM photon antibunching 
SM photon hyperbunching 

Other TM nonclassical effects 



d (n) (l,^) < 
d (,l) (Aal 2 , Aaia) = d^\l,a[ 2 ,a 12 ) < 

d< n >(i, i> ) < o 

d (n) (l,V ) <0 

d< n >(l,wy < -|min((ni),(n 2 )) 

S*\l,W,)<-\\\{C)\-(D)\ 

d (n) (l,n! ±n 2 ) < 

d (n) (A,/2) <0 

d (n) (m,ri 2 ) < 

d (n) (nj -h 2 ) < 

d {n) [h(t),h(t + t)\ < 

d (n) [An(t),An(t + r)] 
= d (n) [l,h(t),h{t + t)} < 

d (n) (l,aia 2 ,a|a 2 ) < 

d (n) (l,aial,a\a 2 ) < 

d (n) (l,ai + a\,a\ + a 2 ) < 

d (n) (l,ai + a 2 ,a\ + 4) < 

d (n) (l,a 1 ,al,al,a 2 ) < 



ED, dA6t 

Ja71- JaToT > 
i TaTH . dABl i 

((AID, dA20t 
( tA2Tl l. JA25T >. jA26b 
l lA3ll . JAM! 

fB4t . dB6t 
l lB7t . dB8t 
fC3l. dC6t 
fC7l . Jcl2t . tC13\ 

GB 



(i) By applying Criterion 3 for F = (1, S1S2, 44)> we obtain 

4 n) = £>((aia 2 ), (o?a|), (nxn 2 )) < 0, (18) 

where ?ii = ajai and 71,2 = 4^2- 

(ii) For = (1, ai4) 0102) one obtains 

df = D((ai4), {a?(4) 2 ), <ftin 2 » ^ 0. (19) 

(iii) For F = (1, ai + a\, a\ + 0,2), Criterion 3 leads to 

df = D«oi + 4), ((a x + al) 2 ), z) ^ 0, (20) 

where z — (hi) + (712) + 2Rc(aia2). 

(iv) For F — (1, a\ + &2, a\ + 4) one has 

df = D(( ai + a 2 ), ({ax + a 2 ) 2 ),z) '< 0, (21) 

where z — (hi) + (^2) + 2Re(di4)- These nonclassicality 
criteria, given by Eqs. ([T8l-(l2"Tli. will be related to the entan- 
glement criteria in subsection lHID 21 



Another example, which is closely related to the Simon 
entanglement criterion ll68ll . as will be shown in subsec- 
tion MI D 1 1 can be obtained from Criterion 3 assuming F = 
(1, 01, a\, a\, 62). Thus, we obtain: 



1 


(Si> 


(4) 


(4) 


(a 2 ) 


(4) 


(a\ai) 




(44) 


(a\a 2 ) 


<&i> 


(a\) 


(ajai) 


(ai4) 


(aia 2 ) 


(« 2 ) 


(&ia 2 ) 


(a\a 2 ) 


(a\a 2 ) 


(al) 


(4) 


(ai4) 


(44) 


<(4) 2 > 


(a\a 2 ) 



(22) 



III. ENTANGLEMENT AND NONCLASSICALITY 
CRITERIA 

Here, we express various two- and three-mode entangle- 
ment inequalities in terms of nonclassicality inequalities de- 
rived from Criterion 3, which are summarized in Table II. 
First, we briefly describe the Shchukin-Vogel entanglement 
criterion, which enables the derivation of various entangle- 
ment inequalities. 
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TABLE II: Entanglement criteria via nonclassicality criteria. 



Reference 



Entanglement criterion 



Equivalent nonclassicality criterion 



Equations 



Duan et al. f661 
Simon [681 

Mancini et al. [69J 

Hillery & Zubairy [67J 

ditto 

ditto 

ditto 

Miranowicz et al. [611 
Other entanglement tests 



d r (Aa 1 ,Aa 2 ) = d r (l, &i, a 3 ) 
d r (l, oi, a[, 02, &l) < 

d r (l, &i + a,2, a\ + &l) < 
d r (l,a 1 a 2 ) < 
d r (l,a?c%) < 

d F (ai,a 2 ) < 
d r (l, 010203) < 
d T (ai, 0,20,3) < 
d r (l,aia l 2 aT)<0 
d T (a k 1 ,a' 2 aT) < 
d r (l, aia 2 ,a\a\) < 
d r (l, oiaj, 0I02) < 
d r (l, &i + a 2 ,&l + a\) < 



<0 dW(Aai,Ao+) = dW(l,ai,ot)<0 g8j-(50j 

d (n) (l, ai, of, 4, a 2 ) + d (n) (l, ai, 4) (52} 
+<f (n) (l, ai, aj, 4) +d (n) (l, ai, al,a 2 ) < 

) d (n \l,a 1 +al,a\+a, 2 ) + 2d < - n ' > (l,a 1 +al) + l <0 (HU, <59) 

d (n) (l,aiaj) < (El, 436} 

d (n) (l,a7(4) n ) < El-El 

d Cn) (ai,aJ) < 01, 07J 

d (n) (l,alo 2 a 3 ) < 03, CU 

d (n) (a|,a 2 a 3 ) < (39) 

d (n) (l,(a|) fe a|a]P) < gl, 

d (n) ((4) fe ,a 2 aD < o gl, gl 

d< n )(l,aiat,ala2) + ({ni+n 2 ) + l)d (n) (l,ai4) <0 4H, 421 
d^{l,aia 2 ,a[al) + (m)(n 2 > + <m + n 2 > d^(l, aia 2 ) < 4H, 4H 

) d (n) (l,ai + 4,4 + a 2 ) + 2d (n) (l,ai +4) < €>2)> 4ED 



A. The Shchukin-Vogel entanglement criterion 

The Criterion 3 of nonclassicality resembles the Shchukin- 
Vogel (SV) criterion ll59l - l6lll for distinguishing states with 
positive partial transposition (PPT) from those with nonposi- 
tive partial transposition (NPT). In analogy to Eqs. (|7]i and 48), 
one can define a matrix M(p) — [Mij(p)] of moments as fol- 
lows: 

M i:j (p) = Tr [(a tn a 12 & ti3 b li ) f (a tjl a? 2 b^ 3 V 1 )p] , (23) 

where the subscripts i and j correspond to multi-indices 
(ix, 12,13,14) and j 2 , j 3 , ji), respectively. Note that, con- 
trary to Eq. 42), the creation and annihilation operators are not 
normally ordered. As discussed in Ref. ||6lll . the matrix M(p) 
of moments for a separable state p is also separable, i.e., 



M(p) 



E 



Pi M A (pt)®M*(pf), 
(24) 

where Pl > 0, J2iPi = M A (p A ) = 

Y^i'ji Mi'j'{p A )\i'){j'\ is expressed in a formal basis 
{\i')} wifhi' = (ii,i a ,0,0) and/ = 01,72,0,0); M B (p*) 
defined analogously. Reference II59I1 proved the following 
criterion: 

Criterion 4 A bipartite quantum state p is NPT if and only if 
M(p r ) is NPT. 

The elements of the matrix of moments, M(p r ) = (p r )}, 
where T denotes partial transposition in some fixed basis, can 



be simply calculated as 

M l0 {p v ) = Tr [(a^a^b^ti^^a^a^b^V 4 )^] 

= Tr [{a^ 1 a^b^V 4 )^ (a^ 1 a j2 b^ 3 b u )p] . (25) 
Let us define 



4(p) 



(/ri fri ) ( /ri fr 2 ) 
(fr N fri) (/r J v/f 2 } 



(/ri fr N ) 
(fr 2 f r N ) 

(fr N f r N ) 



(26) 



in terms of (fiJ rj ) T = ((f%f rj ) r ) = l,...,N). For 
example, if X is an operator acting on two or more modes, and 
we take partial transposition with respect to the first mode, 

X r = (T®id)(X), 

with T the transposition acting on the first mode and id the 
identity operation doing nothing on the remaining modes, re- 
spectively. Then the SV Criterion 4, for brevity referred here 
to as the entanglement criterion, can be formulated as fol- 
lows J6lll: 



Criterion 5 A bipartite state p is NPT if and only if there ex- 
ists F, such that dFp (p) is negative. 

This Criterion 5 can be written more compactly as follows: 

p is PPT O VF : d r p (p) > 0, 
p is NPT & 3F : a%,(p) < 0. 



(27) 
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As for the case of the nonclassicality criteria, the original SV 
criterion actually refers to a set F given by monomials in the 
creation and annihilation operators. This entanglement cri- 
terion can be applied not only to two-mode fields but also 
to multimode fields i6ll l76ll . Note that Criterion 5 does not 
detect PPT-entangled states (which are part, and possibly the 
only members, of the family of the so-called bound entangled 

states) 1 77]. Analogously to the notation of < , we use the 
symbol e < to indicate that a given inequality can be fulfilled 
only for entangled states. 

Here we show that various well-known entanglement in- 
equalities can be derived from the nonclassicality Criterion 3 
including the criteria of Hillery and Zubairy ll67ll . Duan et al. 
ll66ll . Simon l68ll . or Mancini et al. l69Tl . We also derive new 
entanglement criteria and show their relation to the nonclassi- 
cality criterion. 

Other examples of entanglement inequalities, which can be 
easily derived from nonclassicality criteria, include lf78l - [80ll . 
However, for brevity, we do not include them here. 



B. Entanglement and the Cauchy-Schwarz inequality 

The matrix Mp(p) is linear in its state p = J^iPiPi- 
Therefore we have 

M { ;\p) = J2^M i ;\p l )>o (28) 



if M^\pi) > for all pi. Thus, is positive for separa- 

F F 

ble states if it is positive on factorized states. 
Let 



F=(fa...,fa) 



with functions f l = fa fa ■ ■ ■ f iM , where 



(29) 



fa ' _ _ (30) 

I either gj [a,j ) or gj (a! ) if i = kj . 



Here, i is the index of the element fi in F, and index j refers 
to the mode, fa is possibly different from the identity for one 
unique value i = kj, and in that case it is equal to a function 
<7j of either the creation or annihilation operators of mode j, 
but not of both. 

Writing the matrix M^™' in a formal basis {\k)\, one then 

F 



has 



MP = Y.(--flfr-)\k)(l\ 

kl 



(31) 



For factorized states holds 

MP = Y,(-fUn ■■)■■■(■■ flu fm mm 

kl 

= E< : A t iAi :)•••{: fUfkM :)\k)(l\ 

k 

+E( ; A t i^ ■■)■■■(■■ flu fm mm 

= £< : #Jto ■■)■■■(■■ fLJkM :)\k)(l\ 

k 

+ Y,(fii)(fn)---(fL){fiM)\*i)(l\ (32) 

k^tl 

>E(^i>(/*i>-(^><A«>|ft>(i| 

k 

+ T,(fki)(fn)---(fL)(fiM)\*>)(l\ 

k^l 
k 

''£{fii)-{fM){i\)>o. 



c 



The first equality comes from the state being factorized. The 
third equality is due to the fact that the fas are functions of 
either annihilation or creation operators, but not of both, so 
<: fLfli ■■) = (fUn) or (: ftjn :) = </,!&), and that 
for k 7^ I at least one among fa and fa, let us say, e.g., 
fa, is equal to the identity — in particular this implies that its 
expectation value is equal to (fa) = 1. The first inequality is 
due to the fact that (: p k J kl :) = (flJki) or {: flJ kl :} = 
(/fci/fci)' anc ^ to tne Cauchy-Schwarz inequality. 

C. Entanglement criteria equal to nonclassicality criteria 

By applying the nonclassicality Criterion 3, we give a few 
examples of classical inequalities, which can be violated only 
by entangled states. 

1. Hillery •-Zubairy >'s entanglement criteria 



Hillery and Zubairy 167H derived a few entanglement in- 
equalities both for two-mode fields: 



(nin 2 ) °< \(a!al)\ 2 , 
(hi)(h 2 ) <' |(aia 2 }| 2 , 
and three-mode fields 

(nih 2 n 3 ) <* \ (a\a 2 a 3 )\ 2 . 



(33) 
(34) 

(35) 



ki 



These inequalities can be derived from the entanglement Cri- 
terion 5 J59|,|6ll] assuming: F = (1, 0162) to derive Eq. ( 1331 ). 
F = (01, 62) for Eq. ( 1341 ). and F = (1, ^10203) for Eq. 
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On the other hand, Eq. ( f33b can be obtained from the non- 
classicality Criterion 3 assuming F = (1, aia 2 ), which gives 



(n) _ 



1 (aia|) 
(d\d 2 ) (nin 2 ) 



n <5'0. 



(36) 



Analogously, assuming i* 1 = (a%, a 2 ), one gets 



(ni) (aja|) 
(aia 2 ) (n 2 > 



ncl „ 
< U i 



(37) 



which corresponds to Eq. J34t . By choosing a set of three- 
mode operators .F = (1, d{d 2 d 3 ), one readily obtains 



df = 



1 (a}a 2 a 3 ) 
(d\d\d\) (hih 2 hz) 



< 0, 



(38) 



which corresponds to Eq. ( 1351 1. 

By applying Criterion 3 with F = (a\, d 2 a 3 ), we find an- 
other inequality 



(ni) (did 2 d 3 ) 
(aia 2 a 3 )* (n 2 n 3 ) 



< 



0, 



(39) 



which was derived in Ref. 16 lH from the entanglement Crite- 
rion 5. 

Using the Cauchy-Schwarz inequality, Hillery and 
Zubairy rf67ll also found a more general form of inequality 
than the one in Eq. (1331) . which reads as follows: 



((d\raT(4r^) <\(dT(4ry 2 



(40) 



This inequality can be derived from the nonclassicality Crite- 
rion 3 for F = (1, d™(a 2 ) ra ), which leads to 



df = 



i <ar(4r> 



'< 0. (41) 



Alternatively, Eq. (l40b can be derived from the entanglement 
Criterion 5 for F = (1, d™d 2 ). Thus, we see that 



d^(i,d?(dlT) = d r (i,d?d%) c ^o, 



(42) 



where, for clarity, we use the notation d k {F) instead of d k p 
for k = (n) , r. Moreover, we can generalize entanglement 
inequality, given by Eq. (l38l , as follows: 



{n\h\f%) t\({d\fd l 2 d™)\ 



t \k - I p-m\\2 



(43) 



for arbitrary integers k,l,m > 0. This inequality can be 
proved by applying both Criteria 3 and 5: 



dW(i, (&t) k d l 2 &T) = dF{l,a k 1 a l 2 af) 

1 ((a\) k a 2 a?) 
((d\) k a 2 a?)* (n k n l 2 h?) 



ncl „ 



(44) 



where the first mode is partially-transposed. Analogously, 
Eq. ( f39t can be generalized to following entanglement in- 
equality: 



(h k )(h l 2 f%) ^\{d\d l 2 a^)\\ 
which can be shown by applying Criteria 3 and 5: 

d {a l7 a 2 a 3 ) 

(h\) (dXa l 2 a™) 
d k d l 2 af)* (n 2 nr) 



(45) 



dW((a\) k ,a l 2 d?) 



n < 0. (46) 



It is worth remarking that in all the above cases, once the 

< inequalities are found as nonclassicality inequalities, it 
is easy to check that they can be satisfied only by entangled 
states, that is they really are e < inequalities. Indeed, the de- 
terminant condition is the only nontrivial one for establishing 
the positivity of the involved 2x2 matrices. Further, these 
matrices are linear in the state with respect to which the ex- 
pectation values are calculated. Thus, if we prove that the 
matrices are positive for factorized states, then we have that 
they are necessarily positive for a separable state, and so are 
the determinants. For the sake of concreteness and clarity, we 
prove the positivity of the 2x2 matrix of Eq. d37l i for a factor- 
ized state. The positivity of the other matrices for factorized 
states is analogously proved. 

For a factorized state, as a special case of inequalities given 
in Eq. (f32j, we have 



(hi) (d\a\) \ = 
(did 2 ) (h 2 ) I 



> 



(d\di) (d\)(dl) 

(di)(d 2 ) (d\d 2 ) 

$><&!> <aJ><4> 

(di)(a 2 ) (d\)(d 2 ) 



(47) 




<«i> (4) ) > o, 



where the first inequality is due to the Cauchy-Schwarz in- 
equality (X^X) > \(X)\ 2 . 



2. Entanglement criterion ofDuan et al. 

A s harp ened version of the entanglement criterion of Duan 
et al. l66ll can be formulated as follows f59ll : 



(AalAai)(Aa 2 Aa 2 ) °< ((A^Aa^ 2 , (48) 

where Ad.; = d.; — (d^) for i = 1,2. Equation d48l i follows 
from the entanglement Criterion 5 for F = (1, Si, d 2 ) ll59ll or, 
equivalently, for F = (Adi, Ad 2 ). It can also be derived from 
the nonclassicality Criterion 3 for F = (Adi, Ad 2 ). Thus, we 
obtain 



(Aa|Aai) (Ad|Ad 2 ) 
(AdiAd 2 ) (Ad3,Ad 2 ) 



ncl „ 



(49) 
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which corresponds to Eq. (l48l l. Alternatively, by choosing 
F = (1, Si, a 2 ), one obtains 



7 (n) 



1 (61) 
(at) (fix) 
(o 2 ) (aia 2 ) 



(44) 



(50) 



which is equal to Eq. d49l ). Thus, it is seen that this nonclas- 
sicality criterion is equal to the entanglement criterion. More- 
over, the advantage of using polynomials, instead of mono- 
mial, functions of moments in F is apparent. The same con- 
clusion was drawn by comparing Eqs. dA9l ) and dAlOl ) or 
Eqs. (lent and CHI 



(1, a\, a\, a,2, o\). We found that Simon's criterion can be ex- 
pressed as a sum of nonclassicality criteria as follows: 



ai n) {l,ax,a\,a! 



2 , a 2) 



+d^\l,a 1 ,a\,a 2 ) + S n \l,ai,al,a 2 ), (52) 



where d^(l, 61, a\, a\, 0,2) is given by Eq. d22b . More- 
over, dp for F = (1,01,61,02), F = (1,01,62,02) and 

F = (1,01,62) can be obtained from (1221 by analyzing its 
principal minors. Thus, one can prove the entanglement for 
a given nonclassicality by checking the violation of specific 
classical inequalities resulting from the nonclassicality Crite- 
rion 3. 



D. Entanglement criteria via sums of nonclassicality criteria 

Here, we present a few examples of classical inequalities 
derived from the entanglement Criterion 5 and the nonclassi- 
cality Criterion 3 that are apparently not equal. More specifi- 
cally, we have presented in subsection MI Cl examples of clas- 
sical inequalities, which can be derived from the entanglement 
Criterion 5 for a given F\ or, equivalently, from the nonclassi- 
cality Criterion 3 for F 2 equal to a partial transpose of F% . In 
this section, we give examples of entanglement inequalities, 
which cannot be derived from Criterion 3 for F 2 = F[. 

States satisfying Criterion 5 for entanglement must be non- 
classical, as any entangled state is necessary nonclassical in 
the sense of Criterion 1. We will provide specific examples 
that satisfying an entanglement inequality implies satisfying 
one or more nonclassical inequalities. This approach enables 
an analysis of the entanglement for a given nonclassicality. 
The main problem is to express d^, = d r (F) as linear combi- 
nations of some d^- n \F^), i.e.: 



4 



E c ^ (n) (^ (fc) ), 



(51) 



where Cfc > 0. To find such expansions explicitly, we 
apply the following three properties of determinants: (i) 
The Laplace expansion formula along any row (or column): 
detM — J^j (~^Y + ^ Mijfiij, where % is a minor of a ma- 
trix M = (Mij). (ii) Swapping rule: By exchanging any 
two rows (columns) of a determinant, the value of the deter- 
minant is the same of the original determinant but with op- 
posite sign, (iii) Summation rule: If some (or all) the ele- 
ments of a column (row) are sum of two terms, then the de- 
terminant can be given as the sum of two determinants, e.g., 
dot (a + a',b + b'; c, d) = det(a, b; c, d) + det(a', b'; c, d). 



2. Other entanglement criteria 

Now, we present a few entanglement inequalities, which are 
simpler than Simon's criterion, but still correspond to sums of 
nonclassicality inequalities. 

Let us denote the following determinant: 



D(x,y,z,z') 



1 XX* 

x z y 
x y z' 



(53) 



(i) Criterion 5 for F = (1, 61 6 2 , a[d J 2 ) results in 

dF p = D({o 1 al),{o\{dl) 2 ),(h 1 h 2 ) 1 z') °<0, (54) 



where z 1 = ((hi + 1)(^2 + 1))- By using the aforementioned 
properties of determinants, we find that the entanglement cri- 
terion in Eq. (154-b can be given as the following sum of non- 
classicality inequalities resulting from Criterion 3: 



d r p = d( n) (l,6i0 2 ,6l6 2 ) 

+ ((hi) + (h 2 ) + l)d^(l,a 1 ai). 

(ii) Criterion 5 for F — (1, oifiJ,, a\a,2) leads to 



dp 



D((did 2 ), (d 1 a 2 ), z, z) 



/\ ent 

< 



0, 



(55) 



(56) 



where z = (n\h 2 ) + (hi) and z' = (nih 2 ) + (h 2 ). Analo- 
gously to Eq. d55l l. we find that the following sum of the non- 
classicality criteria corresponds to the entanglement criterion 
inEq. d56| >: 



1. Simon's entanglement criterion 

As the first example of such nontrivial relation between the 
nonclassicality and entanglement criteria, let us consider Si- 
mon's entanglement criterion |68]. As shown in Ref. l59ll . 
it can be obtained from Criterion 5 as dp <* for F = 



4 



d (n) (l,dia 2 ,ald 2 ) + (ni)(n 2 ) 

+((ni) + (n 2 ))S n \l,a 1 d 2 ). (57) 



(iii) For F = (1, a\ + a. 2 , a[ + 62), one obtains 

4 = D((oi + 6 2 ), ((Si + a 2 ) 2 >, z, z) c ^ 0, (58) 
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where z = (hi) + (112) + 2Re(ai<22) + 1. Analogously to 
the former cases, we find the relation between the entangle- 
ment criterion in Eq. (l58l and the nonclassicality Criterion 3 
as follows: 

dp = d (n) (l,ai + a 2 ,d\ + d 2 ) 

+2d ( - n \l,a 1 +a 2 ) + l. (59) 

(iv) As a final example, let us consider the entanglement Cri- 
terion 5 for F = (1, ax + a,2, a\ + dJ,)- One obtains 

4 = D{(ai + 4), ((ax + a\) 2 ),z, z') ^ 0, (60) 

where z = (hi) + (712) + 2Re(did2) and z' = z + 2, which 
is related to the nonclassicality Criterion 3 as follows: 

4 = S n) (l,ai +al,a\ +a 2 ) + 2d^(l,ai + aj), (61) 

where d^^l ,cii + a,\,a\ + 0.2) is given by Eq. (120b . and 
d^ (1, a 1 + d\) is given by its principal minor. Equation (l60l 
corresponds to the entanglement criterion of Mancini et al. 
S (see also JH). 

IV. CONCLUSIONS 

We derived classical inequalities for multimode bosonic 
fields, which can only be violated by nonclassical fields, so 
they can serve as a nonclassicality (or quantumness) test. Our 
criteria are based on Vogel's criterion II2411 . which is a gener- 
alization of an alog ous criteria for single-mode fields of Agar- 
wal and Tara 12011 and^more directly, of Shchukin, Richter, 
and Vogel (SRV) i2ll l22ll . The nonclassicality criteria cor- 
respond to analyzing the positivity of matrices of normally 
ordered moments of, e.g., annihilation and creation operators, 
which, by virtue of Sylvester's criterion, correspond to an- 
alyzing the positivity of Glauber-Sudarshan P function. We 
used not only monomial, but also polynomial functions of mo- 
ments. We showed that this approach can enable simpler and 
more intuitive derivation of physically relevant inequalities. 

We demonstrated how the nonclassicality criteria intro- 
duced here easily reduce to the well-known inequalities (see , 
e.g., textbooks Or@|, reviews OrUm, and Refs. |25H4l|| ) 
describing various multimode nonclassical effects, for short 
referred to as the nonclassicality inequalities. Our examples, 
summarized in Tables I and II, include the following: 

(i) Multimode quadrature squeezing i4J] and its generaliza- 
tions, including the sum and difference squeezing defined by 
Hillery Q, and An and Tinh iH S3], as well the princi- 
pal squeezing related to the Schrodinger-Robertson indeter- 
minacy relation ll8"ill as defined by Luks et al. HH). 

(ii) Single-time photon-number correlations of two modes, 
including squeezing of the sum and difference of photon 
numbers (which is also referred to as the photon-number 
sum/difference sub-Poisson photon-number statistics) |@], vi- 
olations of the Cauchy-Schwarz inequality and violations 
of the Muirhead inequality J^HHH], which is a generalization 
of the arithmetic-geometric mean inequality. 



(iii) Two-time photon-number correlations of single modes 
including photon antibunching JH [37| and photon hyper- 
bunching B8il4ltl for stationary and nonstationary fields. 

(iv) Two- and three-mode quantum entanglement inequali- 
ties (e.g., Refs. l66l - l69ll V We have shown that some known en- 
tanglement inequalities (e.g., of Duan et al. ll66ll . and Hillery 
and Zubairy 116711 ) can be derived as nonclassical inequalities. 
Other entanglement inequalities (e.g., of Simon J3) can be 
represented by sums of nonclassicality inequalities. 

Moreover, we developed a general method of expressing in- 
equalities derived from the Shchukin-Vogel entanglement cri- 
terion ll59l l60ll as a sum of inequalities derived from the non- 
classicality criteria. This approach enables a deeper analysis 
of the entanglement for a given nonclassicality. 

We also presented a few inequalities derived from the non- 
classicality and entanglement criteria, which to our knowl- 
edge have not yet been described in the literature. 

It is seen that the nonclassicality criteria based on matrices 
of moments offer an effective way to derive specific inequali- 
ties which might be useful in the verification of nonclassicality 
of particular states generated in experiments. 

It seems that the quantum-information community more or 
less ignores nonclassicality as something closely related to 
quantum entanglement. We hope that this article presents a 
useful approach in the direction of a common treatment of 
both types of phenomena. 
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Appendix A: Unified derivations of criteria for quadrature 
squeezing and its generalizations 

Here and in the following appendices, we present a unified 
derivation of the known criteria for various multimode non- 
classicality phenomena, which are summarized in Table I. 
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1. Multi-mode quadrature squeezing 

The quadrature squeezing of multimode fields can be de- 
fined by a negative value of the normally ordered variance 



nnea by ; 

011 Hi 



(AX ) 2 



< 



with AXfj, = X<f, — (Xf/,), of the multimode quadrature oper- 
ator 



X, 



M 

cf> ^ ^ Cm %m((ftm) 



(A2) 



which is given in terms of single-mode phase-rotated quadra- 
tures 



x 



ra V I 



a m cxp(i(f) m ) + a m exp(-i^ m ). 



(A3) 



It is a straightforward gene ralization of the single-mode 
quadrature squeezing n l25ll . In (IA2t . <fi = (0i, </>m) 
and c m are real parameters. In the analysis of physical sys- 
tems, it is convenient to analyze the annihilation (a m ) and cre- 
ation (a m ) operators corresponding to slowly-varying opera- 
tors. Usually, x TO (0) and i m (7r/2) are interpreted as canoni- 
cal position and momentum operators, although this interpre- 
tation can be applied for any two quadratures of orthogonal 
phases, x m (<f> m ) and x m (4> m + tt/2). 

The normally ordered variance can be directly calculated 
from the P function as follows: 



(A4) 



(: (AJQ,) 2 :) = / d 2 a P(a, a*)[X^(a, a*) - (X^)} 2 , 
where 

X^a.,0.*) = 2J c m (a m e^ 



A I 



(A5) 



m— 1 



and a. = (ai, cxm)- From Eq. (IA4t it is seen that a neg- 
ative value of {: (AX^) 2 :) implies the nonpositivity of the 
P function in some regions of phase space, so the multimode 
quadrature squeezing is a nonclassical effect. This conclusion 
can also be drawn by applying Criterion 3. In fact, by choos- 
ing F = (1, X<p), one obtains 



1 



TO 

>2 

which is the squeezing condition ( lAlb 



TO 



= (: (AX^f :) n ^ 0, (A6) 



By applying the Schrodinger-Robertson indeterminacy rela- 
tion {HI, Luks et al. Ip2ll have given the following necessary 
and sufficient condition for the two-mode principal squeezing 



(Aa\ 2 Aa 12 ) < |<(Aa 12 ) 2 >|, 



(A8) 



(Al) where 



ai2 



&i + a 2 , Adi 



ai2 



(ai2>- 



This condition for principal squeezing can be derived from 
Criterion 3 by choosing F = (Ad\ 2 , AS12), which leads to: 



(Aa\ 2 Aa 12 ) ((Aa 12 ) 2 ) 
{(Aal 2 ) 2 } (Aai 2 Aa 12 ) 



< 



(A9) 



Equivalently, by applying Criterion 3 for F = (l,d| 2 ,ai 2 ) 
one obtains: 



An) 



1 (a{ 2 ) (012) 

(Sl2> (^12) {{ai2) 2 ) 

(4) ((«I 2 ) 2 ) (ni2> 



(A10) 



where 



a\ 2 ai2 



h 2 + 2Rc(a\a 2 ) 



The determinants, given by Eqs. JA91 > and (lAlOb are equal to 
each other and equivalent to Eq. dA8b . This example shows 
that the application of polynomial functions of moments, in- 
stead of monomials, can lead to matrices of moments of lower 
dimension. Thus, the polynomial-based approach can enable 
simpler and more intuitive derivations of physically relevant 
criteria. 



3. Sum squeezing 



According to Hillery 113 311 . a two-mode state exhibits sum 
squeezing in the direction <fi if the variance of 



1 



l " a —id> 1 ~ t ~ t 

(aia 2 e v + a[a 2 e 



satisfies 



where 



<(A^) 2 > < -04), 



(All) 



(A12) 



2. Two-mode principal squeezing 



V 



-(hi +n 2 + l) 



For simplicity, we analyze below the two-mode (M = 2) 
case for c\ = c 2 = 1 and <fr 2 — <fri = tt/2. The two- 
mode principal (quadrature) squeezing can be defined as the 
(^-optimized squeezing defined by Eq. ( 1AU : 



mill <:(A^) 2 :><0. 

4>-<P2— <pl=1T/2 



(A7) 



and h, 



= nt 



for m = 1, 2. As for the case of quadra- 



ture squeezing, d\ and a 2 usually correspond to slowly vary- 
ing operators. Let us denote V x = V(4> = 0) and V y = 
V(4> = tt/2). It is worth mentioning that the operators V x , 
{—Vy) and V z are the generators of the SU(1,1) Lie alge- 
bra. Equation (IA12b can be readily justified by noting that 
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[V x , Vy) = iV z , which implies the Heisenberg uncertainty re- 
lation 

((AV X ) 2 )((AV V ) 2 )>\{V Z ) 2 . 

By analogy with the standard quadrature squeezing, sum 
squeezing occurs when mm{((AV x ) 2 ) , ((AV y ) 2 )} < 
(V z )/2, or more generally if Eq. JA12t is satisfied. We note 
that, in analogy to the principal quadrature squeezing, one can 
define the principal sum squeezing by minimizing ((AV^) 2 ) 
over 6: 



4. Difference squeezing 



As defined by Hillery 03H . a two-mode state exhibits dif- 
ference squeezing in the direction <f> if 



1 



((AW,) ) < -\(W Z )\, 



where 



W A 



{aia\e 1 ^ + a\a,2e l<,> ) 



(A21) 



(A22) 



1 



((AV,) 2 ) < -(V z ) 



(A13) 



Conditions (IA12I ) and (IA13b can be easily derived from Crite- 
rion 3. In fact, by noting that 



1 



{{AV^Y) = (:(AV^:) + -(V Z ), 



(A14) 



the condition for sum squeezing can equivalently be given by 
a negative value of the variance (: (AV^) 2 :). On the other 
hand, by applying Criterion 3 for F = (1, V^), one obtains 



1 ( v 4>) 



(AV,) 2 :) ^ 0, (A15) 



which is equivalent to Eq. (IA121 i. So it is seen that sum 
squeezing is a nonclassical effect — in the sense of Criterion 
1. 

Two-mode sum squeezing can be generalized for any num- 
ber of modes by defining the following A/-mode phase- 
dependent operator 13911 : 




t 



satisfying the commutation relation 



n 



(A16) 



(A17) 



Hereafter j = 1, M and we note that \ {C)\ = (C). Thus, 
multimode sum squeezing along the direction tfi occurs if 



<(AV,) 2 } < 



\(G)\ 



One can find that 



<(AV ) 2 > = (: (AV ) 2 :) + 



1(C) | 



(A18) 



(A19) 



Thus, by applying the nonclassicality Criterion 3 for F 
(1, V^), we obtain the sum squeezing condition 



(: (AV,) 2 -)=df ^0, 
which is equivalent to condition in Eq. ( IA1 8b - 



(A20) 



and W z = — n 2 ). The principal difference squeezing 

can be defined as: 



1 



<(AW^><-KW,)|, 



(A23) 



in analogy to the principal quadrature squeezing and the prin- 
cipal sum squeezing. Contrary to the Vi operators for sum 
squeezing, operators W x = W(<f> = 0), W y = W((f> = n/2) 
and W z are generators of the SU(2) Lie algebra. The uncer- 
tainty relation ((AW x ) 2 )((AW y ) 2 ) > (1/4) | (W z ) | 2 , justi- 
fies defining difference squeezing by Eq. ( IA21I ). One can find 
that 

((AW,) 2 ) = (: (AW^) 2 :) + ±((nx) + <n 2 )). (A24) 
By recalling Criterion 3 for F = (1, W,), it is seen that 



df = (: (AW,) 2 :) 0, 



(A25) 



in analogy to Eq. JA15t . And the condition for sum squeezing, 
given by Eq. dA2U , can be formulated as: 



Oft < — 7T mm \ n il- 

F 2 i=l,2 V ' 



(A26) 



So, states exhibiting difference squeezing are nonclassical. 
But also states satisfying 

\\{ni) - (n 2 )| < {(AW,) 2 ) < i«n x > + <n 2 » (A27) 

are nonclassical although not exhibiting difference squeezing. 
The first inequality in Eq. ( |A27| > corresponds to condition op- 
posite to squeezing condition given by Eq. ( IA2U . 

Criterion 3 can also be applied to the multimode general- 
ization of difference squeezing, which can be defined via the 
operator l40ll : 



K M 



(A28) 



fc=l m=K+l 



for any K < M. For simplicity, hereafter, we skip the limits 
of multiplication in Y[u and [] . The commutation relation 



(A29) 
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where 



1. Sub-Poisson photon-number correlations 



C = + fik) J] nm ~Y[n k + n m ), (A30) 



k m 



justifies the choice of the following condition for multimode 
difference squeezing along the direction </> j40ll : 



((AW^) < 



2 ^ „ \(C)\ 



We find that 



((AW ) 2 } = (: (AW ) 2 :) + ^. 



(A31) 



(A32) 



where 



A I 



D = Y[(l + h k ) Y[ n m +Y[n k + n m ) - 2 JJ ft,. 

ft; m k m j — 1 

(A33) 

By applying Criterion 3 for F = (1, VV^), we obtain the fol- 
lowing condition for multimode difference squeezing: 



The squeezing of the sum (h + = h\ + n 2 ) or difference 
= h\ — n 2 ) of photon numbers occurs if 



(: (An ± ) 2 :) < 0, 



(Bl) 



which can be interpreted as the photon-number 
sum/difference sub-Poisson statistics, respectively 
These are nonclassical effects, as can be seen by analyzing 
the P function: 

(: (An±) 2 :)= J d 2 a P(a, a*)[(\ ai | 2 ± |a 2 | 2 ) - (h ± )]\ 

(B2) 

where a = (ai,a2). Thus, photon-number squeezing im- 
plies the nonpositivity of the P function. The same conclusion 
can also be drawn by applying Criterion 3 for F± = (1, n±), 
which leads to 



1 



(n±) 
(n±> (: :) 



= (: (An±) 2 :) ^ 0. (B3) 



W = (:(AW 4 ) 2 :)<±(\(C)\-(D)), (A34) 



2. Agarwal's nonclassicality criterion 



which corresponds to the original condition, given by 
Eq. ( IA3U . For states exhibiting difference squeezing, the 
right-hand side of Eq. ( IA34l i is negative. In fact, if (C) > 
then 

C - I) = -2 JJ n k m (1 + n m ) - 11 h m ) < 0, (A35) 

k \ m m / 

otherwise 

^-X> = -2m(l + T» fc )-JJn*JIJn m <0. (A36) 

\ k k / m 

It is seen that the difference squeezing condition is stronger 
than the nonclassicality condition d$ <' 0. This means that 

F 

states satisfying inequalities 

'\(C)\-{D)) < <: (AW ) 2 :} < (A37) 



are nonclassical but not exhibiting difference squeezing. 



Here, we consider an example of the violation of the CSI 
for two modes at the same evolution time. Other examples of 
violations of the CSI for a single mode, but at two different 
evolution times, are discussed in Appendix C in relation to 
photon antibunching and hyperbunching. 

By considering the violation of the following CSI: 



(: h\ :)(: h\ :) > (nin 2 ) 2 , 



(B4) 



Agarwal 13 111 introduced the following nonclassicality param- 
eter: 



y/(: h\ :)(:nj :) 

h2 = -p. : r 1. 

(nin 2 ) 



(B5) 



Explicitly, the nonclassicality of phenomena described by a 
negative value of /12 is also implied by Criterion 3 for F = 
(n 1; n 2 ), which results in 



(: hi :) (ftin 2 ) 
(ftin 2 > (: h\ :) 



(B6) 



Appendix B: Unified derivations of criteria for one-time 
photon-number correlations 

Various criteria for the existence of nonclassical photon- 
number intermode phenomena in two-mode radiation fields 
have been proposed (see, e.g., Refs. lUrEHlUlH)- Here, 
we give a few examples of such nonclassical phenomena re- 
vealed by single-time moments. 



3. Lee's nonclassicality criterion 

The Muirhead classical inequality [82;] is a generalization 
of the arithmetic-geometric mean inequality. Lee has formu- 
lated this inequality as follows l34ll 

D 12 = (: hf :) + (: h\ :) - 2(^2} > 0. (B7) 
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The nonclassicality of correlations with a negative value of 
the parameter £> 12 is readily seen by applying Criterion 3 for 

F = (hi — na) = {n~), which yields 



D 12 = (:h 2 _:) ^0. 



(B8) 



For comparison, let us analyze Criterion 3 for F = (1, 
which leads to 



df = (:r\i:)-(h_y>0. 



, 2 cl 



Clearly 



D 12 < d ( - n) ^ 0. 



(B9) 



(BIO) 



Thus, the criterion given by Eq. ( IB9b detects more nonclassi- 
cal states than that based on the D\ 2 parameter. 
Alternatively, a direct application of the relation 

D 12 = J d 2 aP(a,a*)(\ ai \ 2 - \a 2 \ 2 ) 2 n <$' (Bll) 

also implies the nonpositivity of the P function in some re- 
gions of phase space. 



For stationary fields [i.e., those satisfying G' 2 ' (t, t + r) = 
G^(t) so g {2) {t 7 t + t) = 5 (2) (V)], Eq. (Olreduces to the 
standard definition of photon antibunching 0151]: 



. 9 ( 2 )(r)>^(0) 



(2)/ 



(C4) 



Photon antibunching, defined by Eq. (IC3t , is a nonclassi- 
cal effect as it corresponds to the violation of the Cauchy- 
Schwarz inequality: 



G« (t, t)G™ (t + r,t + T)> [G« (t, t + r)\ 



(C5) 



As shown in Ref. 112411 . this property follows from Criterion 3 
based on the generalized definition of space-time P function, 
given by ©. In fact, by assuming F = (n(t), nit+r)), which 
leads to 



( h 2 (t)° ) (°n(t)n(t + t)°) 

(°h(t)n(t + t)°) (°h 2 (t + T)° ) 

GW(t,f) G( 2 )(M + t) 

G^(t,t + T) G^(t + T,t + T) 



2. Photon hyperbunching 



'< 0. (C6) 



Appendix C: Unified derivations of criteria for two-time 
photon-number correiations 

Here, we consider the two-time single-mode photon- 
number nonclassical correlations on examples of photon an- 
tibunching and photon hyperbunching. 



1. Photon antibunching 



The photon antibunching (J, |5|, [3, |g, 14311 of a stationary or 
nonstationary single-mode field can be defined via the two- 
time second-order intensity correlation function given by 

G^{t,t + r) = {°h{t)h{t + r)°) 

= {a"< (t)^ (t + r)a(t + r)a(t)) (Cl) 

or its normalized intensity correlation functions defined as 



9 {2 \t,t + r) 



G (2) (i, t + t) 



y/GW(t,t)GW(t + T,t + T)- 



(C2) 



where ° ° denotes the time order and normal order of field op- 
erators. Photon antibunching occurs if (t, t) is a strict lo- 
cal minimum at r = for g( 2 ' (t, t + r) considered as a func- 
tion of r (see, e.g., Refs. flEl): 

g ( ®(t,t + T)>gW{t,t). (C3) 

Photon bunching occurs if g^ (t, t + r) decreases, while pho- 
ton unbunching appears if g^ 2 '(t, t + r) is locally constant. 



Photon hyperbunching 14 111 , also referred to as photon an- 
tibunching effect ||38t|. can be defined as: 



gW{t,t + T) >gW(t,t), 
given in terms of the correlation coefficient JH 



g {2 Xt,t + r) = 



G (2) (t,t + T) 



G (2 \t,t)G (2 \t + T,t + T) 



(C7) 



(C8) 



7(2), 



where the covariance G (t, t + r) is given by 

G (2 \t,t + r) =G {2) (t,t + r) -G (1) (i)G (1) (i + r), (C9) 

and G (1) (i) = (h(t)) = (ctf(t)a(i)) is the light intensity. 
It is worth noting that, for stationary fields, the definitions 
given by Eqs. ( IC3t and iCli are equivalent and equivalent to 
definitions of photon antibunching based on other normalized 
correlation functions, e.g., 



9 \t,t + T) [ G (l) (t) ]2 ■ 



(CIO) 



However for nonstationary fields, these definitions corre- 
spond in general to different photon antibunching effects J37l 

Analogously to Eq. JC3I >. the photon hyperbunching, de- 
fined by Eq. (IC7b . can occur for nonclassical fields violating 
the Cauchy-Schwarz inequality: 



G {2 \t,t)G {2 \t + T,t + T) > \G m \t,t + t)Y . (Cll) 
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Again, the nonclassicality of this effect can be shown by 
applying Criterion 3 for the space-time P function, given 
by (O, assuming F = (An(t), An(t + r)), where An(t) = 
n(t) — (n(t)). Thus, one obtains 



G (2) (M) 



G (2) (M + t) 



G i2 \t,t + r) G (2) (t + r,t + r) 



'< 0, (C12) 



which is equivalent to Eq. dC7) . Alternatively, by choosing 
F = (1, h(t), h(t + r)), one finds 



d< n) = 



1 (n(t)) <n(* + r)> 

(n(t)) (°n 2 (t)°) (°n(i)n(i + r)°) 

(n(f + r)) (°h(t)h(t + r)°) (°h 2 (t + T)°) 

(C13) 



which is equal to the determinant given by Eq. (IC 1 2b - By 
comparing Eqs. ( IC12b and (IC13l l, analogously to Eqs. iA9\ 
and ( lAlOb . it is seen the advantage of using polynomial, in- 
stead of monomial, functions of moments in F. 



Finally, it is worth noting that the single-mode sub- 
Poisson photon-number statistics, defined by the condition 
(: (An) 2 :) < 0, although also referred to as photon an- 
tibunching, is an effect different from those defined by 
Eqs. (IC3b and dC7b . as shown by examples in Ref. 113511 . 
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